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ASYMPTOTIC BEHAVIOUR OF TOTAL GENERALISED 

VARIATION 


KONSTANTINOS PAPAFITSOROS AND TUOMO VALKONEN 


Abstract. The recently introduced second order total generalised variation 
functional TGV^ ^ has been a successful regulariser for image processing pur¬ 
poses. Its definition involves two positive parameters a and (3 whose values 
determine the amount and the quality of the regularisation. In this paper we 
report on the behaviour of TGV^ ^ in the cases where the parameters a, (3 
as well as their ratio ^/a becomes very large or very small. Among others, 
we prove that for sufficiently symmetric two dimensional data and large ratio 
/3 /q;, TGV| ^ regularisation coincides with total variation (TV) regularisation. 
Total Variation, Total Generalised Variation, Regularisation Parameters, As¬ 
ymptotic Behaviour of Regularisers. 


1. Introduction 

Parameterisation of variational image processing models has not yet been solved 
to full satisfaction. Towards the better understanding of such models, we study 
the behaviour of their solutions as the parameters change. Within the constraints 
of these proceedings, we concentrate in particular on the asymptotic behaviour of 
total generalised variation jBKPlQj . 

In the variational image reconstruction approach, one typically tries to recover an 
improved version u of a corrupted image / as a solution of a minimisation problem 
of the type 

(1.1) min ^(/, Trt) + T(rt), 

u 

where T is a linear operator that models the type of corruption. Here the term 
^{f^Tu) ensures the fidelity of the reconstruction to the initial data. The term 
T(i 4), the regulariser^ imposes extra regularity on u and it is responsible for the 
overall quality of the reconstruction. The two terms are balanced by one or more 
parameters within T. A typical example is T(?i) = aTV(ix), i.e., the total variation 
of u weighted by a positive parameter a jCL97[ IROF92j . While total variation 
regularisation leads to image reconstructions with sharp edges, it also promotes 
piecewise constant structures leading to the staireasing effeet. The second order to¬ 
tal generalised variation TGV^ q, [BKPlOj resolves that issue by optimally balancing 
first and second order information in the image data. The TGV^ q, functional reads 

TGVU (m) = min a\\Du - w\\m + P\\£w\\m, 

wGBD(O) 

where || • \\m is the Radon norm, BD(f]) is the space of functions of bounded 
deformation in the domain f is the symmetrised gradient and [3 > 0. 

Since the values of a and P determine the amount and the quality of the recon¬ 
struction, it is important to understand their role in the regularisation process. In 
this paper we study the asymptotic behaviour of TGV^ q, regularised solutions for 
the extremal cases, i.e., for large and small values of a, (3 and their ratio P/a. For 
simplicity we focus on the case where ^{f^Tu) = ||/ ~ but in most cases, 

our results can be extended to more general fidelities. 
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Summary of our results: In Section 3.1 we show that as long as at least one 
of the parameters o, P is going to zero then the solutions converges to the 

data /. In one dimension we obtain a stronger result, showing in addition that for 
small values of [3 the solutions are continuous. In Section [3^ we focus on the case 
when the ratio P/a\s large, proving that in this regime TGV^q, is equivalent to TV 
modulo “an affine correction”. In Section |3.3| we show that by setting the values 
of a and /3 large enough we obtain the linear regression of the data as a solution. 
In Section [3^ we exploit the result of Section [33] and we show that for sufficiently 
symmetric data and large TGV^ q, is equal to aTV. Our paper is furnished 
with some numerical experiments in Section [33] which verify our analytical results. 


2. Preliminaries and Notation 

In this section we briefly review the basic theory of functions of bounded varia¬ 
tion, properties of TV and TGV^ q, and we also fix our notation. 

Let O be an open, bounded domain in A function u G L^{Q) is a function of 
hounded variation {u G BV(f2)) if its distributional derivative Du is represented by 
an M^-valued finite Radon measure. The total variation of u is defined as TV {u) = 
IIDrIItvt, where ||T||m denotes the Radon norm of an M^-valued distribution T in 

(2.1) lirilx := sup{(r,v) : ||i;||oo < l} , 

and it is equal to the total variation \Du\{Cl) of the measure Du when u G BV(r^). 
The measure Du can be decomposed into the absolutely continuous and singular 
part with respect to the Lebesgue measure Du = D^u + D^u = VuC^ + D^u^ 
where Vu is the Radon-Nikodym derivative D^u/CP. The space BV(f)) is a Banach 
space endowed with the norm ||r||bv(0) = II'^IIl^o) + We refer the reader 

to [AFPOPj for a complete account on the functions of bounded variation. 

Analogously we define the space of functions of hounded deformation BD(r^) as 
the set of all the (O; M^) functions whose symmetrised distributional derivative £u 
is represented by an ^-valued finite Radon measure |TS8Qj . Not at ion-wise one 
can readily check that ||fR||;vi = \£u\{ft). The space BV(f2) is strictly contained in 
BD(f2) for d > 1 while BD(f]) = BV(f2) for one dimensional domains Q. We are not 
going to need much of the theory of BD functions apart from the so-called Soholev- 
Korn inequality. The latter says that if Ct has a Lipschitz boundary then there 
exists a constant Cbd > 0 that depends only on ft such that for every w G BD(f2) 
there exists an element r^ ^ Kerf such that 

(2.2) \\w - r^llLi(o) < CbbW^'^^Wm- 

Here the kernel of £ consists of all the functions of the form r{x) = Ax + 6, where 
6 G and A G is a skew symmetric matrix. 

The second order total generalised variation TGV^^^(r) of a function u G L^{Q) 
is defined as |HKPinirRKVr^lRVn| 

(2.3) TGV^ (m) = min a\\Du - w\\m + P\\£w\\m, 

for (T, /d > 0. The above definition is usually referred to as the differentiation cascade 
definition of TGV^ q,, see |BKPlQj for the original formulation. It can be shown 
that TGV^ Q, is a seminorm and together with || • ||li( 0) form a norm equivalent to 
II ' I|bv(0) |BVllj . i.e., there exist constants 0 < c < C that depend only on ft such 
that for every u with TGV^ q,(r) < oo 

(2.4) c||u||BV(n) < l|w||Li(n) +TGV^_„(u) < C'||M||BV(n). 














ASYMPTOTIC BEHAVIOUR OF TOTAL GENERALISED VARIATION 


3 


Notice that the optimal w in ( |2.3[ ) is not unique in general. In fact re is a solution 
of an problem (not strictly convex). Indeed since + 

we have: 


(2.5) 


re G argmin a\\Du — w\\m 

wGBD(O) 

we argmin / \\/u — w\dx-\ -||fR||x. 

weBD(O) Jq 


In the following sections, we will take specific advantage of the fact that w solves 
(2.5), a problem which can be seen as an analogous one to L^-TV minimisation. 
Let us finally mention that properties of TGV^ q, regularisation have been studied 


in the one dimensional case and when ^(/, Tr) = -|| / — ' 

|BKV18[[PBT^[P^ . 


Ilp(O) 


for p = 1 or 2, in 


3. Asymptotic behaviour 


3.1. /3 ^ 0 while a is fixed and a ^ 0 while /3 is fixed. In this section we 
study the limiting behaviour of regularisation for small values of a, (3. We 

first prove that by fixing a oi ^ and sending or a to zero respectively, then the 
regularised TGV^ q, solution converges to the data /. For simplicity we work on the 
L^-TGV^ Q, denoising problem, i.e., T = Id, but the next result can be extended 
in the more general case e.g. when the fidelity term reads ^||/ — Tr||^, with p > 1 
and T being a bounded, linear operator T : 1/^(0) ^ 1/^(0). For convenience we 
set 

(3.1) {up,a,wp^a)= argmin ^\\f - uWI^/q) + a\\Du - w\\m + I3\\£w\\m- 

ueBV{Q) ^ 

wGBD(O) 


Proposition 1. Let fl C open and bounded and f G fl BY (ft). Then 

(i) Fixing a > 0 we have that \\f — ug,a\\‘L 2 (^Q^ ^0 as (3 ^ 0. 

(a) Fixing P > 0 we have that \\f — ug,a\\‘L‘ 2 (Q) "^0 as a ^ 0. 


Proof, (i) Let e > 0 and {p^} 5 >o be a standard family of mollifiers, i.e., ps{x) = 
6 ~^p{x/6)^ where p G C^(M^), and set fs := ps * /. Because (R/ 3 ,a 5 '^/3,a) is an 
optimal pair in (3.1) by setting u = fs and w = Yfs we have the following estimates, 
for some constant C > 0 


2 11/ - Ul3,a\\l^Q) ^2^\f- «/3,a|lL2(n) + a\\Dup,a “ + P\\£Wl3,a\\M 

< l\\f-fs\\h(n)+mi^fs)\\M 

< l\\f-fs\\Un)+mps*Df\\M 

<l\\f-Mhia)+p]\\Df\\M. 


We set d small enough such that \\f — f 6 \\‘^L 2 (^Q^ < e/2. By choosing < 6e/2\\Df\\Mi 
the result follows. 

{ii) The proof is very similar to the (i) case, by setting u = fs and re = 0, instead. 

□ 


Remark: Of course in both {i)-{ii) cases of Proposition]^ we also get 
\\Dujs^o, — Wjs^o,\\M ^ 0 as well as \\£wjs^o,\\M ^ 0 as ^ 0 or o ^ 0. 
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Another interesting behaviour occurs when ^ 0. In |Vall4j . it is proved that 
for an arbitrary dimension and a fixed a > 0 we have 

||a 4 ^ O 5 as /3 ^ 0. 

However it turns out that in dimension one we are able to prove something stronger, 
provided the data are bounded: 


Proposition 2. Let = (a, 6) C M, / g fl BV(f)) and o > 0. Then there 

exists a threshold /d* > 0 sueh that for every ^ < jS"" we have that 

WD^u^^o^Wm = 0 and W/3^a = 

In partieular this means that for P < 


(3.2) 


1 


Ui3^a= argmin -1|/- M||L 2 (n) +/3||D u\\m- 


weBV(o) 


Proof From the optimality conditions derived in |PB15j . we have that 
solve (3.1) if and only if there exists a dual variable v G i^o(^) 

v" = f - (Cf), -v' G aSgn{Dup^a - (C'a), V G l3Sgn{Dw/3^a) (Cib), 

where for a finite Radon measure jx we define 


Sgn(/i) 


d 


e |Ai|) : ||t;||oo < 1, = p, |mI 


— a.e.| 


Note also that there exists a constant C depending only on ft such that the following 
interpolation inequality holds |Eva TOl Section 5.10, ex. 9] 

\\Dvh 2 ^n) < for all e 


(3.3) 


Observe first that (denoting this dual function v by P/3,a) 

(3.4) ||^'^/ 3 ,a||L 2 ( 0 ) ^ 0 as ^ ^ 0. 

Indeed, from Proposition and condition (Cf) we have that \\D‘^Vi3^a\\L‘^(n) 0 

while from condition (C^) we have that ||p/3,a||oo ^ 0 and thus \\v^^a\\L^{n) ^ 0 as 
P ^ 0. Then we just apply the estimate ( |3.3| ). 

From the fact that we are in dimension one and from (2.4) we have for a generic 
constant C 


||'^/3,a||L=«(0) < C'||'^/3,a||BV(0) < C'( ||'^/3,a ||l 2(0) + (R/3,a)) 

< C{\\f - R/3,a||L2(Q) + ||/||l2(Q) + TGy‘p^^{u^^a)) 

<C{\\f\\mn).TGylM)) 


which in combination with condition (Cf) and the fact that / G I/^(f]) implies 
that \\D‘^ViB^a\\oo < A/". Thus from the Arzela-Ascoli theorem we get the existence 
of a sequence ^ 0 and a continuous function v such that P/3^,a ^ ^ uniformly. 
We immediately deduce using (3.4) that v^^a 0 uniformly as /d ^ 0. But then 
condition {Ca) implies that there must exist a Po such that for every P < Pq we have 
Dujs^a = ^s mcasurcs, i.e., = 0 and since otherwise 

there would exist a point x^^^a ^ with for a sequence 

{Pn)neN converging to 0, a contradiction. □ 


Remark: We believe that the above proof sets the basis for an analogue proof in 
higher dimensions even though admittedly this is a hard task. That would require 
an interpolation inequality for p, divp and div^p analogous to (3.3), as well as a 
proof that the TGV^ q, regularised solution remains bounded, for bounded data /. 
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3.2. Large ratio /3/a. Recall from (2.5) that the optimal re is a solution to a 
IItvi type of problem. This motivates us to study some particular properties 
of the general form of such a problem: 

(3.5) min || 5 (-M;||ii(nRd)+A||£: m;||a<, 5 e A > 0. 

wGBD(S2) 


The next theorem states that if the parameter A is larger than a certain threshold 
(depending only on O) then a solution w of ( |3.5| ) will belong to Kerf. This is 
analogous to the L^-TV problem |CEQ5[ IDAGQQj . where there for large enough 
value of the parameter A, the solution is constant, i.e., belongs to the kernel of TV. 


Proposition 3. Let ft C be an open, bounded set with Lipsehitz boundary, 
g G and Cbb the eonstant that appears in the Sobolev-Korn inequality 

\2.2\ . Then if \> Cbd ci'^d wx is a solution of (3.5) with parameter X, then 


(3.6) 


Wx 


: ms{g) := argmin H^f 

wEKerS 




Proof Since wx is a solution of (3.5), it is easy to check that if r^^ is the element of 


Kerf that corresponds to wx in the Sobolev-Korn inequality then, Wx -='Wx 
solves the following problem: 


' Wx 


(3.7) min IK^r - wx) - + \\\£w\\m- 

weBD(O) ^ ^ 

Indeed, we have for an arbitrary w G BD(f2) 

11(5- 


Wq ~ '^a||li( 0 ;R^) + ^II^'^aIIm < 11^ ~ ('^ + ^WA)llLi( 0 ;Rrf) + A||f (u; + r^(;^)||x, 
with the latter being true since 

Wg - ^a||li( 0 ;R<^) + A||fn;A||M < \\9 - ^IkqojR^) + A||fn;||x, Vn; G BD(f]). 
Since Wx solves ( |3.7[ ), setting Gx := g — wx we have that 

||Ga - kkA||Li(0;Rrf) + A||^^a||a/( < H^A ||li (0;Rrf): 
and using the Sobolev-Korn inequality ||ITa||li( 0 ;R^) ^ C^BDlIf kFA||A 4 we have 


(3.8) 


||Ga — kFA||Li(0;Rrf) 


Gi 


BD 


■||^A||Li(0;Rrf) ^ ||Ga||li(0;R^)- 


A simple application of the triangle inequality in (3.8) yields that if A > Gbd, 
then we must have Wa = 0, i.e., wx = r^x from which (3.6) straightforwardly 
follows. □ 


The notation ms{g) can be interpreted as the median of g with respect to Kerf. 
If d = 1, then this is nothing else than the usual median since in that case Kerf 
consists of all the constant functions. The following corollary follows immediately 
from (2.5) and Proposition It says that for large P/a, TGV^ q, is almost 

equivalent to TV up to an “affine correction”. 


Corollary 4. Let Lt be an open, bounded set with Lipsehitz boundary and let 

a,P > sueh that P/a > Gbd- Then for every u G BV(r^) 

TGV^,„(w) = a\\Du - m£{Vu)\\M- 
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3.3. Thresholds for regression. In this section we show that there exist some 
thresholds for a and [3 above which the solution to the regularisation 

problem is the -linear regression of the data /, denoted by 

/* := argmin \\f - 4>\\\-^(^y 

cf) affine 

We are going to need the following proposition proved in [BVllj : 

Proposition 5 Proposition 4.1]). Let Vt he a bounded, open set with 

Lipsehitz boundary. Then for every 1 < p < d/{d — 1), there exists a eonstant 
CBGviP/c^) > 0; that depends only on Tt, p and the ratio (3/a sueh that 

(3.9) \\u — u^\\Lp{n) < Cbgv(/3/c^)TGV^/,,i(r). 

In the next proposition we show the existence of these regression thresholds for 
d = 2 and also for d > 2 under the extra assumption that the norm of the data 
/ controls the norm of the solution for some p G [d, oo]. 

Proposition 6. Let O C be a hounded, open set with Lipsehitz boundary. Sup¬ 
pose that either 
(i) d = 2 and f G BV(f2) or 

(a) d > 2, / G L^{Q) n BV(f2) and there exists a eonstant C > 0 depending only 
on the domain and p G [d, oo] sueh that ||r||lp(0) ^ C\\f\\];^pm\ for u solution 
to 

then there exist , /3'^ >0 sueh that whenever a > ^ then the solution to 

the L^-TGV^ Q, regularisation problem is equal to /*. 


Proof. Suppose initially that d = 2 and / G BV(f^). Then using the Holder in¬ 
equality along with (3.9) and the fact that any function u G BV(f^) that solves the 
L^-TGV^ Q, problem has a LS norm bounded by a constant C depending only on 
/ and not on a, P 


— 2~ '“llL2(n) 


1„ 

+ :j\\u-u 


*llL 2 (n) + - u){u - u*)dx 


(3.10) 


< 2 II/ “ ^\\h{Q) + C{f)\\u - M*||L2(n) 

< ^11/ - Ah(n) + C(/)Cbgv(/?/«)TGV2/„ ,(«). 


Setting = C(/)Cbgv(1) and = o* we have that if a > and /3 > (3^ we 
can further estimate 

bf-rWlnn) < bf-u\\h^a)+C{f)CBGv{l)TGVl./„.y{u) 


2 

<l\\f 

1 


Ii2(n) + 


< 2ll/-«lli^(n)+TGV^.,„.(«) 

The proof goes through for the case (ii) as well, where the only difference is that 
Holder inequality in (3.10) gives two terms \\u — u'^\\lp(^q^\\u — u'^\\Lp*(n) and ||/ — 
u\\Lp{n)\\u — where p* = pl{p — 1) and oo* := 1. These terms can be 

further bounded using inequality ( |3.9| ) (note that p* < d/(d — 1)) and the fact that 

II'^IIlp(o) < C\\f\\LP(n)- □ 
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More explicit regression thresholds have been given in |PB15j both for general 
and specific one dimensional data /. Let us point out that the condition ||r||lp(o) ^ 
C\\f\\Lp{n) and in particular ||r||oo < C'H/Hoo (which can be derived easily for TV 
regularisation with C = 1), as natural as it may seems, it cannot be shown easily. 
However, if proved, it will also have positive implications as far as the inclusion of 
the jump set of the solution to the jump set of the data is concerned, see |Vall4| . 

3.4. Equivalence to TV for large ratio /3/a and sufficiently symmetric 
data. In Corollary we obtained a more precise characterisation of TCV^ ^, for 
large values of the ratio ^/a. In this section we show that at least for symmetric 
enough data /, TGV^ q, regularisation is actually equivalent to oTV regularisation. 
For the sake of the simplicity of the analysis we assume here that ^4 is a two 
dimensional domain, i.e., C We will also need some symmetry for for the 
time being let be a square centered at the origin. We shall prove the following 
theorem. 


Theorem 7. Suppose that ft is a bounded, open square, eentred at the origin 
and let f G BV(0) satisfy the following symmetry properties: 

(i) f is symmetrie with respeet to both axes, i.e., 

f{xi,X 2 )=f{-xi,X 2 ), f{xi,X 2 ) = f{xi,-X 2 ), for a.e. {xi,X 2 )eft. 

(a) f is invariant under it/ 2 rotations, i.e., f{ 0 ^j 2 x) = f{x), where 07^/2 denotes 
eountereloekwise rotation by i: j2 degrees. 

Then if /3/a > Cbb, the problems 

min - [ If — uTdxTGVj) ^(u) and min - [ If — uTdxaTY(u) 
for p > 1 are equivalent. 

Remark 8. The proof of Theorem is essentially based on the fact that the 
symmetry of the data / is inherited to the solution u and thus to Vr. In that case 
we can show that ms{Vu) = 0 something that shows the equivalence of TGV^^^ and 
oTV. Other symmetric domains, e.g. circles, rectangles, together with appropriate 
symmetry conditions for / can also guarantee that Vr has the desired symmetry 
properties as well. The same holds for any fidelities 4>(/, Tr) that ensure that the 
symmetry of / is passed to u. 


Let us also mention that abusing the notation a bit, by mg{\/u) = 0 we mean 
that zero is a solution of the problem (3.6) with g = Vr. 


of Theorem\^ Since 13/a > Cbd, from Corollary we have that the TGV^ q, reg¬ 
ularisation problem is equivalent to 


(3.11) 


min - / If — uTdx cxllDu — 

neBV(O) pj^' 




Thus it suffices to show that mgiVu) = 0. Since / satisfies the symmetry properties 
(i)-(ii), from the rotational invariance of TGV^^^ |BKPlQj we have that the same 
conditions hold for the TGV^^^ regularised solution u. This also means that Vu = 
(ViR, V 2 R) has the following properties for almost all x = {xi,X 2 ) G ft: 


(3.12) Vir(ti,T2) = Vir(ti, -T2), \/2u{xi,X2) = \/2u{-Xi,X2)j 

(3.13) Vu{x) = -Vr(-x), Vir( 07 ,/ 2 t) = V 2 u{x). 

Recalling that 

ms{Vu) = argmin ||Vr - 

wGKer£^ 
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and that ms{\/u) has the form Ax + 6 it is easy to check, see the following lemma, 
that ms{S/u) = 0. □ 

Lemma 9. Let Q be a square eentred at the origin and suppose that ^ = (^ 1 ,^ 2 ) ^ 
1/^(0; M^) satisfies the symmetry properties 

(3.14) g{x) = -g{-x), gi{xi,X 2 ) = g{xi,-X 2 ), ^ 2 (^ 1 , ^ 2 ) =^ 2 (-^i,^ 2 ), 

(3.15) 9 i{O^I2x) = g2{x), 
for almost every x = (xi,X 2 ) G ft. Then the minimisation problem 

(3.16) 

admits w = 0 as a solution. 


wGKer£. 


Proof. Recalling that Kerf consists of all the functions of the form r{x) = Ax + b 
with A being a skew symmetric function, we have that the minimisation (3.16) is 
equivalent to 


(3.17) 


min / \g(x) — Ax — b\dx, 
Ab Jq 


with corresponding optimality conditions 
f / 9{x) - A0 ^i2 X - b 

jA\9(^) 


, Ot^i 2 X y dx = 0, with Ot^i 2 = 


A0^I2X - b\' 

Using the equalities ^ 2 (^ 1 ,^ 2 ) = gii—X 2 ^xi) and gi{—X 2 ^xi) = gi{—X 2 ^—xi) = 
—gi{xi,X 2 ) we have that ^4 = 0, 6 = 0 solve ( |3.17 ) if 


/ 

JQ 

j 

jQ 


,x)dx = 0 


Onl29{x) 

Q \ \0-^/29{x)\ 

X29lixl,X2) - Xigi{-X 2 ,Xi) 


I 


X 29 l{xi,X 2 ) -Xig2{xi,X2) 


\9{x)\ 


dx = 0 


dx = 0 


In \/gi{xi,X 2 )^ + gi{-X 2 ,xi)^ 
with last equality being true since —g{—x) 


JQ 


\9i{xi,X2)\ 


0 , 


9{x)- 


□ 


3.5. Numerical experiments. In this section we verify some of our results using 
numerical experiments. In Figure we confirm Theorem There, we apply aTV 
and TGV^ Q, denoising with fidelity, to a characteristic function of a disk centred 
at the middle of the domain. Figure [^a) and away from it. Figure [^e). Notice that 
the symmetry properties of Theorem are satisfied for the first case. There, we 
observe that by choosing the ratio fi/a large enough, TGV^^^ and aTV solutions 
coincide. Figures [^b) and [^c) . However, they do not coincide for small ratio 
Figure l^d), see also the middle row slices in Figure j^i). In this case TGV^ q, 
produces a piecewise smooth result in comparison to the piecewise constant one of 
aTV. Note that when the symmetry is broken, aTV and TGV^ q, solutions do not 
coincide even for large ratio fi/a, Figures [3g),[3h) andj^j). 

Figure [^depicts another example of an image that satisfies the symmetry prop¬ 
erties of Theorem 7 The aTV solution coincides with the TGV| q, one for large 
ratio fi/a^ Figures ]2T) and[^c), but not for small ratio. Figure |^d). 


Finally in Figure 


we solve the L^-TGV^ q, regularisation problem in a noisy 
image. We observe That for very small values of /3 or a, essentially we have no 
regularisation at all, see Figures |^c) and |^d) respectively, verifying Proposition 
In Figure l^e), we choose a large ratio fi/a^ obtaining a TV-like result which 
is nevertheless quite different than the aTY result. Figure [^f), having staircasing 
only inside the ellipse. This is due to the “affine” correction predicted by Gorollary 
1^ see also the corresponding diagonal slices in Figure |^i). Figure |^g) depicts 
a typical TGV solution with no staircasing while in Figure |^h) we set a and fi 
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(a) Symmetric data (b) TV solution, (c) TGV solution, (d) TGV solution, 
a = 10 a = 10, ^ = 10® a = 10, ^ = 200 



(e) Non-symmetric (f) TV solution, 
data q; = 10 


(g) TGV solution, (h) TGV solution, 
a = 10, /5 = 10® a = 10, /3 = 10® 


— / 

- TGV: a = 10, = 10® - 

- TGV : Q = 10, /3 = 200 



(i) Gorresponding middle row slices for 
symmetric data 



(j) Gorresponding middle row slices for non- 
symmetric data 


Figure 1. Illustration of the two dimensional aTY and TGV^ q, 
equivalence for symmetric data when P/a is large enough. Notice 
that the equivalence does not hold once the symmetry is broken. 



(a) Original image (b) TV solution, (c) TGV solution, (d) TGV solution, 
a = l a = 1,(3 = 100 a = l,(3 = 2 


Figure 2. Illustration of the two dimensional (tTV and 
equivalence for symmetric data when P/ais large enough. 


large enough and we obtain the linear regression of the data, as expected from 
Proposition!^ 
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KONSTANTINOS PAPAFITSOROS AND TUOMO VALKONEN 




(e) TGV solution, 
a = 0.1, (3 = 100 



(f) TV solution, 
a = 0.1 



(g) TGV solution, 
a = 0.1, (3 = 
0.15 


(h) TGV solution, 
a = 100, h = 
1000 



Figure 3. L^-TGV^ q, denoising for extremal values of a and (3. 
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